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Abstract 

We investigate the recent conjecture that the chiral phase transition in non- 
compact lattice QED is driven by monopole condensation. The comparison of an- 
alytic and numerical results shows that we have a quantitative understanding of 
monopoles in both the quenched and dynamical cases. We can rule out monopole 
condensation. 
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1 Introduction 



In a series of papers [1-10] we have investigated strongly coupled QED, both on the lattice 
in the non-compact formulation and in the continuum using Schwinger-Dyson equations. 
The strong coupling region is of interest because of the existence of a second order chiral 
phase transition. This implies a continuum limit, and the interesting question is whether 
or not it describes an interacting theory. Our calculations of the renormalized charge, e#, 
and fermion mass, tur, demonstrated that whenever ttir goes to zero in lattice units (i.e. 
the ultraviolet cut-off is removed) then e R goes to zero. This suggests that the theory 
is non-interacting in the continuum limit in accordance with the general belief that non- 
asymptotically free theories are trivial. It is encouraging that the two approaches, namely 
lattice and Schwinger-Dyson, agree with each other. Further support of this picture comes 
from other authors UTTI. 



However, this picture has been queried by Hands, Kocic, Kogut and collaborators 
[12-20] who investigated the behaviour of magnetic monopoles near the phase transition. 
Using a new monopole 'order parameter' first introduced by Hands and Wensley [^(J , they 
conclude that monopoles condense in the chirally broken phase. The occurrence of this 
proposed second order monopole phase transition is important because such a transition 
would imply the existence of monopoles in the lattice model's continuum limit. This 
would cast doubt on conclusions about continuum physics drawn from lattice calculations, 
as monopoles are presumably absent in QED. Furthermore, dual superconductivity and 
charge confinement are to be expected whenever monopoles condense |[L5|| . This is in 



conflict with our picture [[10] where we have found free electrons and massless photons in 
the broken phase. 

In the strong coupling limit (3 — > 0, or the limit when we have a large number of 
flavours, the action is dominated by the fermion determinant, which is a compact object 
in the sense that it only depends on the compactified link variables e lAlJ -, being the 
gauge field. It was found that the compact U(l) Wilson action has a first order phase 
transition |21] which is driven by monopole condensation |22| . Hence, it is conceivable 



that monopoles play a role in the non-compact well. Indeed, simulations with 



very large numbers of fermions p3| suggest that the phase transition becomes first order. 



In this paper we shall investigate the relevance of monopoles for the phase transition. 
We solve the quenched case analytically |23| and look at the dynamical fermion case 



numerically. We find no evidence of monopole condensation. 
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2 Lattice monopoles 



The action for non-compact lattice QED with dynamical staggered fermions is 

5 = f E W + E^lPw + ^kW, (2-i) 

x,fi<u x,y 



with 



F„ v {x) = A^A u {x) - A^x), (2.2) 
D*,v = ^E(- 1 ) a ' 1+ ^ + ^ 1 [e^ W W-e-^ (?/) W], (2.3) 

where A M is the lattice forward derivative and /3 = 1/e 2 . 



To define monopoles |2T|, ^(J we decompose i 7 ^ into an integer valued string field N pu 
and a compact field / M „ which lies in the range (— 7r, 7r]: 

= 2ttN pv + f^. (2.4) 

The Bianchi identity tells us that F pv summed over any closed surface always gives zero. 
This does not apply to the N pu and f pv fields separately. This allows the common defini- 
tion of a monopole current 



M^x) = ^-e pupa A u f pcT (x + ft) = -\e pupa A u N pa (x + /2), (2.5) 

47T I 



which lives on the elementary cubes of the lattice (or equivalently on the links of the dual 
lattice). Equation (|2.5| ) shows that a string can only end on a monopole or antimonopole. 
Each component of M p can take the values 0, ±1,±2. The current M p is conserved: 
A p M p [x) = 0, where A M is the lattice backward derivative. 

Later on we shall be interested in the monopole susceptibility. This is defined by [^] 

1 4vr 2 
Xm = ££<Wz)«0)> = — £<iV(:r)iV(0)>, (2.6) 

" x ^ x 

where we have made use of the fact that ^2 x F pu (x) = 0. In the infinite volume limit 
further manipulations lead to the equivalent form 

4vr 2 

x m = -—Y.(x 2M ^) M ,m- (2-7) 



The physical interpretation of eq. ( |2.7| ) is that it measures the fluctuations of the total 
dipole moment, whereas eq. ( |2.6| ) can be regarded as the residue of the photon pole in 



the compact photon propagator. We have used eq. (|2.6|) for our measurements to avoid 
ambiguities in defining x 2 on a finite lattice. If the monopoles condense we would expect 
X m to diverge. 

At a second order phase transition fluctuations of physical quantities become large and 
should show up in other quantities as well. We have looked at fluctuations in monopole 
charge. The corresponding susceptibility is 





Figure 1: The susceptibility (2.8) for an ideal Bose-Einstein gas, 
at various temperatures above and below the condensation temper- 
ature T c . 
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It has long been known that when Bose condensation occurs the tendency of bosons to 
occupy the same state causes long range density correlations, which in turn leads to a 
divergence of \ q at small momentum. For an example see 'Distance Correlations and 
Bose-Einstein Condensation' |]26| which finds that in a Bose gas with a condensation 
temperature T c the density-density correlation drops off exponentially if T > T c (giving 
a finite x at p = 0), drops off like r~ 2 at T c (giving a x diverging like and drops off 
like r _1 below T c (giving a x diverging like p~ 2 ). These behaviours are illustrated for the 
case of the ideal Bose gas in fig. 1. We see that condensation gives a spectacular signal. 

One subtlety is that for periodic boundary conditions Xq{p) vanishes for p = on finite 
lattices and so one must extrapolate to zero momentum. (Note that the limits have to be 
taken in the correct order: volume goes to infinity before momentum goes to zero.) 

Other quantities we shall also look at are the monopole density 

P = 777EI M m(^)I (2-9) 

(V: lattice volume) and the string density [] 

ff = IT7 E 1^)1- (2-10) 

We would expect that these quantities would show non-analytic behaviour (though not 
necessarily a divergence) around a phase transition. 



3 Analytic results 



The authors of |20], [17j claim that even in quenched non-compact QED there is a phase 
transition at which the magnetic monopoles condense. Because in the quenched case the 
action is Gaussian, we can derive analytic formulae for most quantities p4j . The simplest 
quantity is the string density. The probability distribution for a single F^ v field is a 
Gaussian: 

y(F) = Tr-^e"^ 2 . (3.1) 



1 The normalization is chosen to agree with ref. [Eo| 
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The distribution is completely determined because we know the width (F 2 ) = 1/(2/3). 
This gives (cf. eq. (fnOQ ) 

= ^(\N(F)\) 



/oo 
dF^{F)\N{F)\ (3.2) 



3 <^ 
2 



oo 

oo 



' ^ erfc ((2n + 1)tt/?5 



^ n=0 

On a finite lattice with volume V the only change is that the width of the Gaussian is 
reduced to (F 2 ) = (V — l)/(2V/3). This means that a on a finite lattice can be found by 
making the replacement (3 — > (3V/(V — 1) in eq. ( |3.2|) . One sees that finite size effects are 
already negligible on rather small lattices (which holds for all the other quantities in this 
section as well). 

The next quantity of interest is the monopole density. To find the monopole density 
we need the probability distribution for the six F fields on the faces of a cube. The 
outwardly directed 'plaquettes' are labelled using the dice convention, namely that F n 
and F 7 _ n are on opposite faces. We find 

(Fi, . . . , F 6 ) = n-h^a - 6)§(a + b)^5{F 1 + ■ ■ ■ + F 6 ) 

(3.3) 

xexp{-(3a(F 2 + ■ ■ ■ + F 2 ) - 2pb(F 1 F 6 + F 2 F 5 + F 3 F 4 )}. 

This is the most general Gaussian form consistent with cubic symmetry and the Bianchi 
identity. (Because F\ + • • ■ + F 6 = we can add an arbitrary multiple of (F\ + ■ ■ • + F 6 ) 2 
to the exponent without changing the distribution \I/ at all. This freedom has been used 
to eliminate terms of the form FiF 2 etc.) The parameters a and b are fixed by the known 
expectation values 



{Ft) 



d 4 k 1 - ci 1 



(3 J (27r) 4 4- Cl -c 2 -c 3 -C4 2/3' 
(2l 

2 r d A k -ci + cic 2 47 - 1 



IFF\ 1 f d ^ k - l + c i + C2 ~ c lC2 7 
{FlF2) ~2pJ (2tt)M - Cl - c 2 - c 3 - c 4 -"27' (3 ' 4) 

(FFe) = 



/9 J (2tt) 4 4- Ci -C 2 -C 3 -C4 2/3 ' 

where c M = cosfc^. (On a finite lattice the integrals are to be replaced by sums over allowed 
non-zero momenta.) On an infinite lattice 

7 = 0.215 563.... (3.5) 
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(Note that the value of 7 is calculated from the correlation between two F fields, which is a 
gauge invariant quantity. Therefore \I/ is gauge invariant.) To give the correct expectation 
values, a and b must take the values 



1 1 + 7 

a = ' - = 0.826 049.. 

127(1-27) 

1 1-57 

b =—— -4- =-0.052 879. 

127(1-27) 



(3.6) 



The monopole density p(/3) is 
p(P) =(\M ll \) 

/oo 
dF 1 ■ ■ ■ dF 6 *(Fi, .., F 6 ) \Ni + ■ ■ ■ + N 6 \. 
-00 



(3.7) 



The fact that |M M | is bounded is enough to show that all derivatives of the monopole 
density p(f3) are finite at all (3 values. If p{(3) is expanded as a series of the form 

00 

p(/3)=/3 | E a n(&-/ 3 ) n ( 3 - 8 ) 

n=0 



Figure 2: The monopole density p and the a to p ratio as a func- 
tion of (3 for quenched QED. The symbols represent the data from 
ref. [20] together with our data, while the curves show the analytic 
results on an infinite lattice. 
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about an arbitrary point /5 , then the bound on |M M | leads to bounds on the a n : 

< a „ < \ ^±|l A -»-. (3.9) 
3 n!m + 1)! 4™ 



These bounds are strong enough to show that the series in eq. (|3.8| ) is convergent with a 
radius of convergence of (at least) j3o- A convergent series expansion rules out the existence 
of any essential singularities in p. There is certainly no sign of a phase transition in p((3). 
A similar proof holds for correlation functions involving a finite number of /'s and iV's. 



In fig. 2 our formulae are checked against the Monte Carlo data ||20||. The agreement 



Figure 3: The same as fig. 2 for a larger range of (3 on a logarithmic 
scale. The data points at the smaller and larger (3 values are our 
measurements. 
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is excellent. In fig. 3 we show the monopole density over a larger (3 range. This enables 
us to show the asymptotic limits of our formulae. At /3 — 0, p{(3) goes to 7/15, while for 
large (3 the density drops off like 6 erfc(7r/?2) ; i.e. approximately exponentially. Also in 
this extended range we find good agreement with the data. The ratio a(f3)/ p(f3) goes to 
1/4 at large (3. This must be so because at large f3 the only topological excitations that 
occur are isolated plaquettes with |F| > tt, and these are surrounded by a monopole loop 
of length 4. At small (3 the ratio diverges like /3~z (cf. fig. 3). In the interesting region 
around [3 = 0.24, a((3)/p(/3) has only grown to ~ 0.34, indicating that monopoles and 
antimonopoles are on a short leash (i.e. are tightly bound). If we look at a timeslice the 
average length of string joining each monopole-antimonopole pair is simply 4<r//9, and is 
only about 1.4 lattice units. 

We now turn to the discussion of the susceptibilities. To calculate these we need 
to know the general two-'plaquette' distribution ty(Fi,Fj) and so the compact photon 
propagator (ftfj). If 

(F?) = (Ff) = -L and <F«*» = || (3.10) 



with pij given by the non-compact photon propagator fl27| , 

(F a ,(0)F, u (x)} = J-/ 7^4 1 " e 



ik-x 



X 



2(3 J (2tt) 4 4 - ci - c 2 - c 3 - c 4 
(6 ait {l - e"^)(l - e*-) - Ml - e-^Xl - e ifc «) 

-5^(1 - e~ ifc «)(l - e <fc ") + 6 Pv (l - - e ifc ")) (3-11) 

then 

m, Fj) = §(1 - PIT 1 * eM-j^iF? + Ff) + ^F iFj }. (3.12) 

" Pij Pij 

With the help of the Fourier series 

oo (_i)k 

f = -2j2—r L MkF) (3.13) 

fc=i K 

we derive the series 

,,,, , ,fi ,fi (-1)"+" . ^(nm Vij \ ( n 2 +m 2 \ 

for the compact propagator. Knowing (fifj) for any pair of plaquettes we can calculate 
Xm from its definition eq. ( |2.6| ). 



In fig. 4 we show the theoretical curve for \ m as a function of (3. At finite (3 there is 
no divergence: XmiP) is a correlation function involving /'s (or A's) for which the proof 
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Figure 4: The monopole susceptibility Xm for quenched QED as a 
function of (3 on an infinite lattice. 

of analyticity applies. At (3 — 0, Xm has the value |7r 2 because the / fields are completely 
uncorrelated and evenly distributed in (— 7T, 7r]. It is striking that after /3 ~ 0.1 the curve 
drops exponentially in /3 because the monopole density is dropping so quickly. 

If we had monopole condensation we would expect that the monopole correlation 
function (x 4 = t) 

C(t)=J2(M,(x,x 4 )M lM (0)) (3.15) 

X 

would exhibit long-range order pB| . C(t) can be evaluated using the general form ( |3.14| ) 
for (fifj)- In fig. 5 we show C(t) on an infinite lattice at the (3 value {(3 = 0.244) where 
Kocic et al. [171 place a phase transition. We find that the correlation function drops by 
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approximately three orders of magnitude between t — 1 and 2. 

The Fourier transform of C(t) gives the charge susceptibility Xqip) ( see ec L- (EH)) when 
the momentum p is directed along a lattice axis. In fig. 6 we show Xq(p) as a function of 
p for three different f3 values. Because we can work on an infinite volume we can go all 
the way to p = 0. We never see the enhanced long wavelength fluctuations characteristic 
of Bose-Einstein condensation (see fig. 1), on the contrary long wavelength fluctuations 
are strongly suppressed. The charge susceptibility does not diverge but vanishes like p 2 
for all (3 values. As usual the finite size effects are very small, and the curve measured on 
a finite lattice is already very close to that seen on an infinite lattice. 



Figure 5: The analytic result for the monopole correlation function 
C(t) for quenched QED at (3 — 0.244 on an infinite lattice. The 
solid (open) symbols represent positive (negative) values of C(t). 
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Figure 6: The susceptibility Xqip) as a function of momentum on 
an infinite lattice for three values of (3 for quenched QED. Kocic 
et al. place a phase transition at /3 — 0.244. The momentum was 
chosen along a lattice axis. 

Many other quantities are also calculable. Particularly simple to calculate are the 
compact Wilson loops, which again are analytic and give a potential between charges 
which is Coulombic at all f3 > 0, inconsistent with confinement. 

The authors of ref. |15j claim that strongly coupled dynamical non-compact QED 
undergoes monopole condensation and confines electric charge, the only evidence they 
present being the existence of a percolation threshold. Their case is greatly weakened by 
the fact that the quenched theory has a similar percolation threshold |17j (in fact, the 
evidence for the divergence of the Hands and Wensley cluster susceptibility is strongest 
for the quenched theory), but as we have seen there is no monopole condensation or 
confinement. 
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4 The dynamical case 



We shall now investigate what happens if dynamical fermions are included. For the action, 
which corresponds to four fermion flavours, see eqs. (2.1-2.3). 

If, as seems likely from the last section, the monopole properties are determined by 
very short distance fluctuations of the electromagnetic fields, we could expect that these 



Figure 7: The monopole density as a function of the plaquette 
energy. The data are for four flavours of dynamical staggered 
fermions. The solid symbols represent our data, while the open 
symbols represent data from ref. [19]. The curve is the analytic 
result for the quenched case on an infinite lattice. 
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properties are determined by the plaquette energy values, because the plaquette energy 
is a good measure of the fluctuation strength. 

Therefore we have plotted the monopole density against P = ^"^^{F^J) in fig. 7 
using data from refs. [|T^, [K| and our results given in Table 1. We also show the analytic 
curve calculated in the quenched case. We find surprisingly good agreement between the 
data and the analytic result, indicating that the inclusion of dynamical fermions does 
not change our previous conclusions. The data comes from a wide range of bare masses, 



Figure 8: The cr to p ratio as a function of the plaquette energy. 
The data are for four flavours of dynamical staggered fermions as 
given in Table 1. The triangle is m = 0.04, the square is m = 0.02 
and the circle is m = 0.01. The curve is the analytic result for the 
quenched case on an infinite lattice. 
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m = 0.005 — 0.16, and plotting against P has brought them all on a universal curve. (Note 
that the quenched case can also be viewed as the m — > oo limit.) 

In fig. 8 we show the ratio cr/p. The measured values are about 2% higher than the 
quenched calculation. Plotting against P has again brought measurements at different 
masses onto the same curve. The fact that a is a little different in the four flavour case 
than in the quenched case shows that the single F distribution deviates slightly from a 
Gaussian form. 



Figure 9: The monopole susceptibility Xm as a function of the pla- 
quette energy. The data symbols are the same as in the previous 
figure. The curve is the analytic result for the quenched case on an 
infinite lattice. 
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The monopole susceptibility \ m is plotted against P in fig. 9. Again, we compare the 
data with the analytic result. As has been remarked before ||21||, this quantity is hard to 



measure accurately because of large cancellations between positive and negative charges. 
Within the errors we find agreement with the analytic quenched result, and do not see 
any divergence of \ m . (The 'phase transition' reported in is at P ~ 1.025.) 

In fig. 10 we show Xq{p) as measured on our four-flavour configurations compared 
with curves calculated in the quenched case at the same value of P. Measurement and 
calculation are in excellent agreement. Even at small (3 (the upper curve) Xg vanishes at 
small momenta, showing that the monopoles have not condensed. Consider for example 
X q measured at the momentum p = 2tt/L. This measures charge density fluctuations 
with a wavelength equal to the size of our lattice. Since a wave with wavelength equal 
to the lattice size is positive in one half of the lattice and negative in the other we are 
measuring the difference between the monopole concentrations in the two halves of our 
lattice. If the monopoles have truly condensed the tendency of bosons to occupy the same 
quantum state would cause large differences between the number found in the right-hand 
and left-hand sides of our lattice. This happens for a condensed boson gas (see fig. 1), but 
for monopoles we see that long wavelength fluctuations are strongly suppressed, showing 
that condensation has not taken place. 



Figure 10: The susceptibility Xqip) as a function of momentum 
for three values of j3. The data are for four flavours of dynamical 
staggered fermions at m = 0.04. The curves are analytic results for 
the quenched case on an infinite lattice. 
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The surprising agreement between all the theoretical curves (with no adjustable pa- 
rameters) and the Monte Carlo data confirms our hypothesis that the monopole properties 
are determined by short-distance fluctuations. 





m 


P 


P 


a 


Xm 


0.17 


0.04 
0.02 


1.2832(6) 
1.2669(6) 


0.2087(2) 
0.2041(2) 


0.07621(9) 
0.07395(9) 


2.729(107) 
2.816(89) 


0.18 


0.04 
0.02 


1.2022(6) 
1.1881(6) 


0.1865(3) 
0.1826(3) 


0.06533(13) 
0.06355(12) 


2.371(150) 
2.813(146) 


0.19 


0.04 
0.02 
0.01 


1.1343(5) 
1.1194(5) 
1.1106(9) 


0.1673(2) 
0.1626(2) 
0.1603(3) 


0.05662(7) 
0.05459(7) 
0.05355(11) 


2.320(82) 
2.198(88) 
2.050(122) 


0.20 


0.04 
0.02 
0.01 


1.0739(3) 
1.0617(4) 
1.0548(6) 


0.1493(2) 
0.1457(2) 
0.1437(2) 


0.04901(7) 
0.04745(6) 
0.04666(10) 


2.110(85) 
1.969(77) 
2.095(121) 


0.21 


0.04 
0.02 
0.01 


1.0234(4) 
1.0133(4) 
1.0068(4) 


0.1342(2) 
0.1314(2) 
0.1292(2) 


0.04289(6) 
0.04176(6) 
0.04091(9) 


1.890(63) 
1.806(64) 
1.899(89) 


0.22 


0.04 
0.02 


0.9779(3) 
0.9692(3) 


0.1205(2) 
0.1179(2) 


0.03766(6) 
0.03670(5) 


1.595(56) 
1.476(51) 



Table V. The monopole density p, the string density a and the 
monopole susceptibility Xm on a 12 4 lattice with four flavours of 
dynamical staggered fermions. Also given are the plaquette energy 
values P = i En<v( F L)- 



5 Discussion 



In this work we have investigated monopoles in quenched and dynamical non-compact 
lattice QED. In the quenched case we have derived analytic formulae which we have 
checked against numerical data. Here we can prove that there are no singularities and 
divergences in the quantities we have looked at. The same formulae describe the dynamical 
case when quantities are plotted against the plaquette energy, which measures the strength 
of the electromagnetic field. Thus we have arrived at a quantitative understanding of 
monopoles in both the quenched and dynamical case. 

The similarity of the percolation threshold in the quenched, Nf = 2 and Nf = 4 
cases ||24j| , occurring at almost the same monopole density and with almost identical 
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critical exponents, lends further support to the idea that the monopole behaviour is the 
same in the quenched and dynamical case. This picture may, however, change when Nf 
is so large that the phase transition becomes first order p3 . 



We have also looked at the distribution of monopoles and antimonopoles in a time 
slice. We saw that 60% of all monopoles have an antimonopole on adjacent sites and 
only 10% a monopole. So actually what we see looks more like a gas of dipoles than a 
condensate of monopoles. (One can also see this from the low cr/p value; see figs. 2,3, 8.) 



In dynamical QED we have already checked [10] that the potential is Coulombic and 



the photon does not acquire a mass. This is also inconsistent with confinement at low (3. 

It is not surprising that when the monopole density becomes large (p ~ 0.15) perco- 
lation takes place. However, percolation is not necessarily connected with condensation 
or with any other field-theoretic or thermodynamic property of the theory. Indeed, it is 
rather easy to find examples where the percolation threshold and the "authentic" phase 
transition are at different couplings. One example is the Ising model of higher dimension, 
where the percolation threshold lies at higher f3 than the phase transition f28| . Also, if one 



looks at the 3d Ising model at (3 = and non-zero magnetic field, h, one finds percolation 
thresholds at certain values of h, because if (3 = the spins are randomly distributed 
with the concentration controlled by the magnetic field. It is well known that randomly 
distributed sites on a cubic lattice first percolate when the concentration has reached 
~ 32% PUL So for strong negative h the up-spins do not percolate, while the down-spins 



do. There is a percolation threshold at about 32% (h ~ —0.38) and then both spins 
percolate. At a concentration of ~ 68% (h « +0.38) there is a second threshold above 
which only the up-spins percolate. Despite the occurrence of these percolation thresholds 
there are certainly no phase transitions at non-zero h in the Ising model. 



The cluster susceptibility of Hands and Wensley [^(J is 

_ / (S n =4 9n n ) ~ n max \ /p- i\ 

where n is the number of dual sites in a cluster linked together by monopole world lines, 
g n is the number of clusters of size n and n max is the size of the largest cluster. This 
susceptibility has long been used to find percolation thresholds (it is essentially the Sip) 
of H] or the x f of 0). However it can lead to misleading results if it is used to locate 
phase transitions due to the following defects. First of all, Xc is not a Green's function, so 
a divergence of Xc does not imply an infinite correlation length and so does not indicate 
a second order phase transition. Furthermore, Xc counts monopoles and antimonopoles 
with the same sign, whereas physically they should contribute with opposite sign. 

Moreover, in order to compute Xc one must treat monopoles differently depending on 
whether or not they belong to the same cluster. No physical operator can do this because 
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by Bose symmetry all physical operations treat indistinguishable particles identically If Xc 
is not a physically realizable quantity the fermions cannot couple to it and the divergence 
of Xc cannot cause a chiral phase transition. 

It is worth noticing that in dynamical non-compact QED Xc diverges in places where 
there is no phase transition. The chiral phase transition takes place only at m = 0, while 
at finite m quantities such as (xx) are smooth functions of j3. However, Xc diverges not 
only at m = but for all m including m = oo, (see ref. fl8|, [H 



In conclusion, the papers [12-20] on monopoles in non-compact QED do not prove 
that monopoles are relevant in the continuum limit of the lattice theory, and so do not 
invalidate the picture of the chiral phase transition presented in refs. [1-11] . 
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